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Abstract Interferometric calibration always yields non unique solutions. It is there- 
fore essential to remove these ambiguities before the solutions could be used in any 
further modeling of the sky, the instrument or propagation effects such as the iono- 
sphere. We present a method for LOFAR calibration which does not yield a unitary 
ambiguity, especially under ionospheric distortions. We also present exact ambigui- 
ties we get in our solutions, in closed form. Casting this as an optimization problem, 
we also present conditions for this approach to work. The proposed method enables 
us to use the solutions obtained via calibration for further modeling of instrumental 
and propagation effects. We provide extensive simulation results on the performance 
of our method. Moreover, we also give cases where due to degeneracy, this method 
fails to perform as expected and in such cases, we suggest exploiting diversity in time, 
space and frequency. 

Keywords Instrumentation: interferometers - Techniques: interferometric - 
Methods: analytical 



1 Introduction 

Self calibration is essential for radio interferometers such as LOFAR0 to obtain high 
quality results under the presence of corruptions. In this paper, we focus most of our 
attention on propagation effects caused by the ionosphere. The ionosphere is an active 
medium which affects electromagnetic radiation passing through it. Naturally, radio 
interferometric observations on the earth are distorted by its effects due to this reason. 
Most significant distortions were first observed in very long baseline interferometry 
(VLB I) due to the fact that extremely long bas elines will see the sky through unrelated 
parts of the ionosphere. We refer the reader to lCotton ( 1995 1. (and references therein) 



for an overview of earliest attempts of VLBI calibration under the influence of the 
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ionosphere. Such attem pts have become simp ler since the introduction of the matrix 
measurement equation ( Hamaker et all 1996 ). which encapsulates the full polarized 
form of radio interferometric data in a concise way. 

Ionospheric distortions can also become significant at very low frequencies, even 
with baselines of moderate length. Therefore, instruments such as LOFAR, MWA, 
PAPER, GMRT, and SKA which can observe at frequencies close to ionospheric cut- 
off will certainly have to cope with such disturbances. Significant work has already 
being done in calibration and correction for ionospheric distortions of i nterferomet- 
ric observations at very low frequency and the reader is referred to e.g., (Ilntema et al 
20091: Ivan der Tol and van der Veenll2007h . This paper is not about tackling such dis- 



tortions, but rather about extraction of ionospher ic (as wel l as ins trumental) infor- 
mation with minimal ambiguities. As shown by lHamakei d2000b . self calibration 
always introduces certain ambiguities of the solutions. At this point, we shall clar- 
ify the difference between a degeneracy and an ambiguity. A degeneracy arises when 
we do not have enough observations to extract complete information from the data. 
For example, consider locating a point in three dimensional space, just by distance 
measurements from a set of known reference points. If we have only two distance 
measurements, the location of the point will have a degeneracy. That is, it could lie 
anywhere on the intersection of the two spheres (with radii equal to the measured dis- 
tances) from the reference points. However, by having three or more properly placed 
reference points, we could eliminate this degeneracy. In contrast, an ambiguity arises 
due to our description of a problem, or due to the choice of a coordinate system to de- 
scribe our parameters. Take the case of finding the square root of a given number: We 
know that we could always have more than one (complex) solution to this problem. 
Note that this cannot be eliminated by taking more measurements. 

In the case of radio interferometry, we have more than enough baselines compared 
to the number of stations (or unknowns) to eliminate any degeneracy. However, the 
ambiguities remain, regardless of the number of baselines or data points we have. 
This is not a major concern in most cases with an unpolarized sky as the ambiguity 
will not affect the end result. On the other hand, to construct models for instrumen- 
tal and propagation effects we need to use solutions obtained by calibration. In such 
situations, it is essential that we use ambiguity free solutions. Traditional radio tele- 
scopes such as the Giant Meterwave Radio Telescope (GMRT) and the Very Long 
Baseline Array (VLBA) eliminate the ambiguities by periodically using a noise in- 
jection source but this if far from practical for a telescope like LOFAR. 

Calibration of a phased array interferometer such as LOFAR has its own chal- 
lenges as well as advantages in comparison with traditional dish based (movable) 
interferometers. Due to not having any movable components, approaching electro- 
magnetic radiation from a celestial source will have a time varying relationship with 
the LOFAR dipoles because of the rotation of the sky. On the other hand, due to the 
same reason (no movable parts), the LOFAR beamshape is simpler to simulate using 
numerical electromagnetic software. In contrast, a steerable dish based telescope can 
track a source across the sky but the beamshape is more complex to simulate due to 
shadowing, spillover etc. Numerical simulation of LOFAR element beamshape has 
given us extensive knowledge of its behavior. Therefore, in this paper, we shall use 
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the a priori knowledge of the LOFAR beamshape to our fullest advantage to restrict 
ambiguities in calibration. 

This paper is organized as follows: In section|2] we give a brief overview of LO- 
FAR calibration and the arising ambiguities. Next, in section|3] we derive expressions 
for the ambiguities in closed form. In section |4] we discuss the necessary conditions 
to obtain reduced ambiguity solutions and cases where these conditions will not hold. 
Finally, we present results based on simulations in section [5]bef ore concluding. 

Notation: Matrices are denoted in bold uppercase letters (e.g. A) and vectors in 
bold lowercase letters (e.g. v). The sets of Complex, Real and Integer numbers are 
denoted as C,R and Z respectively. An estimate (or a possible solution) for a given 
quantity (e.g. ff) is given by 0. The column vector with all ones is denoted by 1 and 
all zeros by 0. The submatrix created by selecting a subset of rows (or columns) from 
matrix A is represented as A.[ a ^ C :d\, where we have selected rows a to b and columns 
c to d. The pseudoinverse of A is denoted by A' while the norm by | All- 



2 Calibration 



Consider the correlations observed on baseline p-q: We can express th e observed vis- 



ibility matrix (ignoring additive noise) for baseline p-q in the standard (IHamaker et al 
1996) Jones formalism as 



Y pq - GpE p F p CF^E^G^ 
where we have N stations and (p, q) are drawn from the set 

5 = {(l,2),(l,3),...,(l,A0,...,(iV-l,A0}. 



(1) 



(2) 



The number of baselines (cardinality of S) is Nb < N(N-l)/2. The matrices G, E and 
F in ([TJ correspond to the instrumental gain, element beamshape and the ionospheric 
Faraday rotation, (qualified by subscript p or q to denote the station) respectively. 
We first make the following assumptions, mainly to simplify our analysis: 

- We observe a dominant, unpolarized, point source, at the phase center (therefore 
no Fourier phase component) and the station beams are pointing at that source 
(therefore station beam gain is unity and E represents the dipole element beam). 

- We assume electronic leakage is small (therefore G is diagonal). Moreover, if 
there is any leakage it is assumed to be part of the dipole E Jones matrix. 

- We assume we have a priori knowledge of the E Jones matrix. Moreover, because 
all stations have parallel dipoles (as in LOFAR), the E Jones matrices are assumed 
to be identical for all stations. 

- Ionospheric delay phase is absorbed into the clock delays in G Jones. 

Note that the assumptions made here are specific for a typical LOFAR observation 
and therefore the results based on these assumptions are mainly applicable to a tele- 
scope that is very similar to LOFAR. 

The unknowns for station p are given in G ; , and F^: 



?, P 
gyp 



cos(a'p) sin(a p ) 
- sin(a' ; ,) cos(ap) 



(3) 
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where g X p,gyp G C and a p e R are the unknown station gain and Faraday rotation 
parameters, respectively. The element beam shape determines the E (= E p — E 9 ) 
Jones matrices, about which we have fairly accurate knowledge from numerical sim- 
ulations. The source coherency (with intensity /) is given as 



/ 
/ 



(4) 



Consider the problem of calibration: We can represent (0 as 

\„ = J P Cjf (5) 

where J p and J q are the Jones matrices that we need to estimate. However, we can 
get solutions for (0 such as J p = J P XJ, where U is an unknown unitary matrix. This 
is perfectly fine if we only need to correct the data using these solutions because the 
unitary ambiguity cancels out in the final product. However, if we need to study the 
parameters in (0, it is hard to disentangle the unitary ambiguity from the solutions. 

Another way of looking at this problem is by counting the number of degrees of 
freedom that is actually present in the problem. If we consider the unknowns g X p,gyp 
(4 real numbers), a p (1 real number), we get 5 degrees of freedom per station. How- 
ever, if we attempt calibration as in (0, we have 8 degrees of freedom, which is 
clearly an overkill. 

Therefore, instead of solving for a system as in (0, we use the original form 
with unknown parameters as given in (0. Still, this does not guarantee us an 
ambiguity free solution as we see later. Hence, we study the possible ambiguities of 
(0 in section0and give the necessary conditions to get a solution in section|4] 



3 Ambiguity Analysis 



Consider the unknowns g X p,gyp in (0. We can rewrite this as g xp = \g xp \exp{j(/> xp ), 
gyp = \g yp \ exp(j<p yp ) where <p xp , <p yp are the unknown (true) phases. In order to study 
the ambiguities we could get for the unknowns in (0, we could rewrite the possible 
solutions as 



(-l) k '"\g X p\x 
exp [j(cf> xp + k xp n + k ap n)j 









(-i)M 



exp (j((f>yp + kypit + k ap n)) 



(6) 



and 



cos(a p + k ap n) sm(a p + k ap n) 
- sin^r, + k ap n) cos(a p + k ap n) 



(7) 



where k xp , k yp , k ap e Z. Note that if we substitute (O and instead of (0 into 0, 
we should still get the same value for \ pq . The necessary conditions for us to rewrite 
(0 as (0 and will be discussed in section^ 
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Consider the product 



E pq = EpFpCFqEq 



(8) 



such that V E 



= GpEpqG^ . Assume E pq to have non zero off diagonal elements 



i.e., Ei2,pq + 0, JZvL, m + 0. The necessary conditions for this is related to the same 
conditions for © and O to be valid solutions to (Q]i and will be given in section|4] 

Then, in order to get the same \ pq in (Q]), by substitution of G p , G q given in (O 
as well as G p , G q given in © to (Q~|), we get the following conditions: 



+ kyp = 0, 2, . . . 

kxp kyq — 0, 2, . . . 

kyp + k xq — 0, 2, . . . 

kyp + fcyg 0, 2, . . . 



(9) 



and 



4>x P - 4>xq + k xp n - k xq n + k ap n - k aq n ■ 
4>x P - 4>yq + k xp n - k yq n + k ap n - k aq n ■ 

4>yp ~ 4>x q + ky p 7T - k xq JT + k op 7T - k aq JT ■ 

(py P - 4> yq + kyp-n - k yq n + k ap n - k aq n - 



•rxp 

<t>xp 
4>yp ' 
Pyp ' 



(pyq 
<Pxq 
byq. 



(10) 



We can write similar condition for all other baselines. 

Accumulating all possible constraints such as (0 for all baselines, we have 

Ak = n 

AN b x2N 2Nxl 4N b xl 



(ID 



where k = [k x \ , k y \ , k X 2, k y 2 ■ ■ ■ , k y pj] T . The vector n has only even integers as entities. 
This is a linear Diophantine equation (because k e Z 2N ). The matrix A^ has ones at 
two locations on each ro w. We can show this has rank of 2N (full column rank) using 
similar techniques as in ( Chen and Petropulul2001 ; Yatawatta et all 2004 ). Therefore, 
we get a unique solution to k (ideally equal to 0). 

Next, we can also accumulate equations in (TTOb for all possible baselines to arrive 

at 



AN b x2N 



4> + k n + k tr n 

y 2Nx\ ZNxl 2Nxl . 



= A n 



(12) 



AN b x2N 2Nxl 



where = \$ xX , <p yX , cf> x2 , <p y 2, <pNy\ T , 4> = [0*1 - <Py\ . 4>x2, 4>y2, <pNy] T , and k a = 
[k a \,k a \,k a 2, k a 2, ■ ■ . , k a m] T . The matrix A p has 1 and -1 at two loc ations on each 
row. Similar as in dChen and Petropulu 2001 ; Yatawatta et al 2004), we can show 
that A p has rank 2N - 1 (null space of rank 1). Therefore, in order to find a solution 
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to (fT2l i. we have to fix one element in 0, say by making <f>^ = <p x \ = 0o- Then, we 
can rewrite (fTSl i as 

A-p[ % 2:2N] [^>[2:2N] + k[2:22V]^ + ka[2:2iV]7Ij (13) 
= Ap [: ,2:2iV]^[2:22V] _ A p[:.l] (~0[1] + 00 + k[i]7T + k a [i]^J 



which gives us the solution 



(14) 



9[2:2N] - 9[2:2N] 
- (k[2:2iV]7r + Ki[2:2N] n ) ■ 

To simplify dl41) . we need to find A^. 2 -2N] A pt- l ~l - ^ e I 310 * that tne sum °f eacn row 
of A p is zero. Therefore, 

A p l = (15) 

A^[: ; 2:22V]1 + A p [ : ,l]l = 

Therefore, we can rewrite ([14} as 

0[2:2iV] = 0[2:2/V] + 1 (~0[1] + 00 + k[i]7T + k ff[1] 7r) (16) 
- (k[2:2ATr + Kr[2:2N]ri) • 

Closer scrutiny of ([TBI reveals that the estimated phase, 4>[2-2N] i s equal to the 
true phase, 0[2:2iV] with a constant unknown phase (common to all stations) and an 
ambiguity of kn, where k e Z. Therefore, the ambiguity of the solutions for the X 
and Y phases for a given station is the same (within n). The difference between the 
X phase and Y phase ambiguities is called the phase zero difference in traditional 
calibration jargon. Thus, we get no phase zero difference (except an ambiguity in n) 
for all the stations. 

Next, we shall investigate the conditions for F p , F ? and F p , F q to satisfy (fTJ to get 
the same Y pq . We arrive at equations similar to 

a p - a q + k ap n - k aq n = a p - a q (17) 

for each baseline. Combining all possible equations (fTTl i for all baselines, we get 



/ ^ 
Ay ^_®^ + ka 7T 

N,,xN V Nxl Nx\ / N b xN 



Ay (18) 

Nxl 



where or = \a\,a>2, ■ ■ ■ ,aN] T , a = [a\, 02, ■ ■ ■ , (XnY and k s = [k a \, k a 2, . . . ,k a ti\ T - 
The matrix A y has 1 and - 1 at two locations on each row. We can once again show 
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that the rank of this matrix is N — 1 . Hence, in order to get a solution, we fix one 
element of S, say S[i] = Si = oq. Then, we can rewrite ([18} as 

A/[:,2:An (a![2:N] + ^a[2:N]^) = ^f[:,2:N]^[2:N] (19) 
-A/t :> i] (-«[!] + Q-o + k5[i]7r) 

to yield a solution as 

S[2:W] = CH[2:N] - A^ 2:JV] A/r : ,i] (-£*[l] + ffo + k s[ i]7r) - k 5 [2 : Arj7r. (20) 

Once again, because the sum of each row in Af is zero, using similar analysis as (Q3), 
we can show that AL 2 -at|A/[:,1] = _ 1- Therefore, we can simplify the solution as 

a { 2:N] = a[2:N] + 1 ( _Q, [1] + a + kff[l]7r) - k s[2: iV]7r. (21) 

Therefore, our estimate ap-.N] is equal to the true values, a^-.N] plus an unknown 
constant which is common to all solutions with an ambiguity of kn, k G Z. Thus, we 
can estimate the differential Faraday rotation within a common unknown rotation and 
an integer ambiguity of n. 



4 Necessary Conditions 

In this section, we give necessary conditions for our approach to work, especially 
with respect to the properties of the element beam pattern or the E Jones matrix. 



4. 1 Non diagonal E pq in © 



First, we study the requirements for 
Jones matrices in <[8J as 



to be non diagonal. We express the E and F 



Ep - E 9 - 



E x e E xt f, 
E y e £ V 



(22) 



cos(a p - a q ) sin(a p - a q ) 
- sin(ar p - a q ) cos(a p - a q ) 



where E x g, E^, E y g, E yi p are the field components in the directions (9, <p in cylindrical 
polar coordinates. The conditions required for the off diagonal terms to be zero then 
becomes: 

cos{a p - a q XE xe E* e + ) (23) 

+ sin^ - a q ){E xe E* y4> - E x4 ,E* B ) = 

cos(a ; , - a q )(E* e E ye + E^E y< j) 

+ sin(Qr p - a q )(-E* e E yi p + E^Eyg) = 0. 

Note also that the off diagonal terms of the product EE H is equal to E x gE* g + E^E*^ 
(and its conjugate). If this term is zero, we get a solution to ( 1231 as a p — a q + 2kn. 
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On the other hand, if ExgE*, - E xt f,E* g = 0, then we get another solution to ( l23l l as 
a p — a q + (2k + l)7r/2. So, if either of these conditions are satisfied by the element 
beam, our method fails to give a unique solution. 

A more general way of looking at (l23T l is as follows: We can represent ( l23l ) in 
matrix form as 

Ec = (24) 



where 



E 



{E x8 E* g + E^E*^) (E x gE^ - E x4 ,E* g ) 
(E x8 E* g + E^E*^* -(ExgE*^ - E x<p E* g Y 

cos(a p - a q ) 
sin(a p - a q ) 



(25) 



c 



We can find a solution for c in (l24l if at least one of the eigenvalues of E is zero (and 
the solution for c is the scaled eigenvector). Therefore, when E has at least one zero 
eigenvalue, it is possible (under certain values of Faraday rotation) that our method 
fails. On the other hand, when none of the eigenvalues are zero, we are guaranteed 
that regardless of the Faraday rotation, our method works. 



4.2 Convex Optimization 

Normally, the solutions for the calibration are found using a non linear optimization 
technique. Therefore, its natural to consider the solution of (Q} as an optimization 
problem. The residual error of ([I] can be given as 

r„ = vec(V pq ) - t pq (G) (26) 

where f pq (0) = vec(G p E p FpCF^E^G^J. The unknown parameter vector is given 
by 9. For each station we have 5 (real) unknowns, i.e., \g xp \, \g yp \, <p xp , <Pyp and a p . 
However, since we can only find differential Faraday rotation (not the absolute ro- 
tation), we keep a for one station fixed. So, the parameter vector can be given as 

= [\g x i\, \ga\, . . . , \g y i |, \g y 2\, 4>xi . <f>x2, • ■ • , (f>yi , (f>y2, ■ ■ ■ , «2, «3, • • -f ■ Hence, for 
stations, the size of 6 is 5A^ - 1 . Finally, we can write the calibration problem as 

= arg min ^ rf }q r pq . (27) 

PA 

The Hessian of this optimization problem will have rank 5A^ - 2 at the solu- 
tion. This is because, we still have an ambiguity in (p xp or <p yp . Therefore, the smallest 
eigenvalue of the Hessian will be zero (or close to zero under noi se). For a well formu- 



lated optimization problem, we need to show that d27i i is convex (IBovd and Vendenber^he 



2004) with respect to 6. This implies that the Hessian should be positive (semi) def- 
inite. In other words, all eigenvalues (barring the smallest eigenvalue which is zero) 
should be positive. Therefore, by studying the second smallest eigenvalue, we can 
study the performance of our calibration, in particular with variation of the beam in 
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different directions of the sky. The more positive this is, the more well formulated our 
problem becomes. 

Under noise, the Cramer-Rao bound is inversely proportional to the second small- 
est eigenvalue that we discussed. Therefore, we can study the same behavior by using 
the Hessian instead of the Cramer-Rao bound in this case. 



4.3 Frechet Derivative and Condition Number 

We have assumed a priori knowledge of E in our analysis. However, there will cer- 
tainly be small variations in the real beamshape from our assumed model. We need 
to show that small changes in E in ([TJ only leads to small changes in V. This ensures 
that the proposed method works, even when the real E will be slightly different from 
our model and moreover, when each station will have a slightly different beamshape. 

We use the Frechet derivative of ([1} to examine the sensitivity of V to small 
changes in E. The Frechet derivative is a generalization of the derivative of scalar 
functions to higher dimensional functional spaces. More details a bout the Freche t 
derivative and sensitivity analysis of matrix functions can be found in lHighaml(l2008l) . 
The dependence of E on Y pq in ([D can be written as 

V(E) = \ pq = GpEFpCFf E H Gf (28) 



Let AE denote a small change to E. Then, as in lHighaml d2008l) . we write the Frechet 
derivative as 

L V (E, AE) = V(E + AE) - V(E) - o(||AE||) (29) 

where 

lim — — = 0. (30) 
IIAEii-o ||AE|| 

Therefore, we get the Frechet derivative as 

L V (E, AE) = GpEFpCF^ AE H Gf + G ;) AEF p CF^ E H Gf . (31) 
The relative condition number of V(E) is given as 

,, v _, ||L V (E,AE)1| ||E|| 

cond(X(E)) = max . (32) 

IIAEIMO ||AE|| ||V(E)|| 

Using norm inequalities ||AB|| < ||A||||B|| and ||A + B|| < ||A|| + ||B|| we get an upper 
bound for the relative condition number as 

cond(\(E)) < ||G /; ||(||F ; ,CF«E H || + ||EF p CFf ||)||G,|| jJjjjjL (33) 

We can use ( 1331 to study the sensitivity of the data to variations in E and hence, the 
sensitivity of our calibration to errors in our knowledge of E. Generally, it is better to 
have a condition number as small as possible, the smallest being equal to 1 . 
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4.4 Diversity 

In this section, we have given several ways of looking at the necessary conditions 
for our approach to work. In practice, there will always be cases where some of 
these conditions are not satisfied. However, even if these conditions are not satisfied, 
we can exploit diversity to make our method work. Diversity can be described as the 
opposite of degeneracy. Indeed, degeneracy arises when we do not have enough inde- 
pendent observations. In contrast, when we do have observations taken independently 
from one another we create diversity. There are several different forms of diversity 
that we can exploit. Implicitly, we used the knowledge of E and moreover, the fact 
that this is not an identity matrix, in our proposed method. In other words, we used 
polarization diversity arising due to X and Y dipoles having different beam patterns. 
Furthermore, even if the aforementioned conditions are not satisfied for one baseline, 
there are other baselines where the ionosphere is different and a solution is possible. 
Secondly, ionosphere will almost always change with time and even for a baseline 
which has degeneracy, with time, things will change to yield a solution. Thirdly, the 
element beam changes with direction in the sky and as the source moves across the 
beam, it will overcome degenerate cases. Moreover, ionospheric properties have well 
known frequency behavior. In situations where we fail to obtain a solution, we can 
extrapolate solutions obtained at different times,baselines, frequencies etc., to get a 
complete solution. 



5 Simulation Results 

We present some results based on numerical simulations in this section. We select 
the observation frequency to be at 120 MHz. The E Jones matrix is derived from 
numerical electromagnetic simulation of a LOFAR high band antenna (HBA) dipole. 
The elements of the E Jones matrix at this frequency is shown in Fig. Q] and Fig. 
12 Both these figures show the beam over the full hemisphere, orthographically pro- 
jected onto the plane. The dipoles are rotated by ji/A radians from the meridian. The 
center of each figure correspond to the zenith. The perimeter of the figure correspond 
to the horizon. Note that in all figures in this section we use the same projection and 
coordinate system (except in Fig. 0(b)). In Fig. [3] we have also shown the power 
beam for both X and Y dipoles. 

Now, let us reconsider the necessary conditions discussed in section [4] In Fig. 
IH we have shown the conditions when ( l23l l could be satisfied. We have shown the 
products ^E^E^ + E^E*^ and |(£ , xe £ , * -£ , A - £ , *)| in Fig.g](a) and (b), respectively. 
The first term goes close to zero along the planes of the dipoles and at the horizon. 
However, the second term only goes to zero at the horizon. 

In Fig. 3](c), we have shown magnitude of the smallest eigenvalue of E (normal- 
ized by ||E||) given in ( 1251 ). This figure tells us something more than the previous two 
figures: which is that apart from the directions along the planes of the dipoles (which 
also includes the zenith), no possible Faraday rotation exists such that ( l23l l is satis- 
fied. In other words, regardless of the Faraday rotation, we should be able to estimate 
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l 2.5 

■ 1.5 



(a) (b) 

Fig. 3 The power beam for X and Y dipoles at 120 MHz. (a) X power beam \E^g\ 2 + \E Ja t\' 2 (b) F power 
beam \E yB \ 2 + \E y4 ,\ 2 



it in all directions except the planes of the dipole and the horizon, because ^ is not 
diagonal. 

So far, we have studied the characteristics of the beam that enables our calibra- 
tion to work. Now, we shall consider the full optimization problem and evaluate its 
performance. We have simulated an array with N - 54 LOFAR stations. We vary the 
direction of the celestial source in the sky (keeping the flux at 1 Jy), and for each 
direction, we do a simulation. The direction is varied at discrete intervals with az- 
imuth range [0, 2n) radians and elevation range [0, n/2] radians. The gains 
are simulated from a uniform distribution in [0.5, 1.5]. The phases <p xp ,(f> V p are simu- 
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-1.6 -1 -0.5 0.5 1 1.5 



(C) 

Fig. 4 Conditions for (8) to be not diagonal, (a) amplitude \(E x gE* g + E^E*^ (b) amplitude \(E x gE*^ - 
ExtpE* s )\ (c) amplitude of smallest eigenvalue of E 

lated from a uniform distribution in [0, 2n) radians. The Faraday rotations a p are also 
simulated from a uniform distribution in [0, 2n) radians. We have also added white 
Gaussian noise to the simulated visibilities, keeping the sig nal to noise ratio at 10. 
We use a stand alone optimization routine (lLourakisll2004l) to solve the calibration 
problem (l27l l. The initial values for Ig^Ugvpl are kept at 1 and the rest of parameters 
are kept at 0. 




(a) (b) 



Fig. 5 (a) The second smallest eigenvalue of the Hessian of )27t . (b) Residual normalized error of 1 1271 . 

In Fig. |5](a), we have shown the second smallest eigenvalue of the Hessian ob- 
tained after calibration. This is positive for most directions in the sky and where it 
goes close to zero, we need to exploit diversity to obtain a complete solution. How- 
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ever, the overall behavior shows that the problem is mostly convex. In Fig. 0(b), we 
have shown the normalized residual error of (l27i . which is small and mostly domi- 
nated by noise, except for some high values at the horizon. 

In Fig [6] (a), we have shown the relative condition number obtained using ( 1331 ). 
The average value of the condition number over all baselines is shown. It is below 
3 for this particular case, indicating that the problem is well conditioned and quite 
robust to error in modeling of E. Next, in Fig. 0(b), we have shown the histogram 
of the errors in estimating (f> xp and <p rp . The true error is small and hardly visible in 
this figure. The dominant feature is the ambiguity in n. We see similar behavior in the 
solutions for a p . 




Fig. 6 (a) Relative condition number )33t . (b) Histogram of the error of the solutions for <f> xp and <j>y p , 
divided by n. 



6 Conclusions 

We have shown that it is indeed possible to calibrate LOFAR with reduced ambigui- 
ties. The analysis was done for an unpolarized point source, and it is straight forward 
to extend this to an unpolarized extended source. The crucial factor in our analysis 
was that the a priori knowledge of the element beam shape. Extensive electromag- 
netic simulations has in fact given us that knowledge. We have also shown under 
which conditions our calibration fails to give an answer. However, in such situations, 
we can exploit diversity to get a complete solution. The extension of this method to 
polarized sources as well as results based on application of the proposed method to 
real LOFAR data will be presented in future work. 
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